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Abstract 

We study MHV amplitude for the 2 — > 5 scattering in the multi-Regge kinematics. The Man- 
delstam cut correction to the BDS amplitude is calculated in the leading logarithmic approxima- 
tion (LLA) and the corresponding remainder function is given to any loop order in a closed integral 
form. We show that the LLA remainder function at two loops for 2 — > 5 amplitude can be writ- 
ten as a sum of two 2 — > 4 remainder functions due to recursive properties of the leading order 
impact factors. We also make some generalizations for the MHV amplitudes with more external 
particles. The results of the present study are in agreement with all leg two loop symbol derived 
by Caron-Huot as shown in a parallel paper of one of the authors with collaborators. 
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1 Introduction 



We further study the high energy behavior of the scattering amplitudes in M = 4 supersymmetric 
Yang Mills (SYM) theory. The recent impressive advances in the field of the Maximally Helicity 
Violating (MHV) amplitudes motivated us to continue the programm of applying the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) jlj approach to the MHV amplitudes, initiated by two of the authors 
in collaboration with A. Sabio Vera in Ref. [2,3. In Ref. [2 the analytic properties of the Bern- 
Dixon-Smirnov (BDS) ansatz [4] were tested and found to contradict the well known high energy 
behavior of the multileg scattering amplitudes starting at two loops for six external particles. A 
similar conclusion was drawn from the strong coupling side by Alday and Maldacena [5]. The 
analytic properties of the BDS amplitude were also studied in Refs. [5] [71 [5]. 

Two of the authors with a collaborator in Ref. 3 derived a closed integral form of the all 
loop correction to the n = 6 BDS formula in the leading logarithmic approximation (LLA) and 
calculated it analytically at two loops. The need of correcting the BDS amplitude by so-called 
remainder function triggered a big activity in this direction, which resulted in introducing novel and 
powerful computation techniques. The most recent one is the use of the symbol successfully applied 
by Goncharov, Spradlin, Vergu and Volovich [5] to calculating the six-gluon two-loop remainder 
function from Wilson loops with light-like edges |T0l [IT] . Their result was shown by two of the 
authors [12] to reproduce the BFKL predictions in the multi-Regge kinematics and was used to 
calculate the next-to-leading impact factor p~3] needed in the BFKL approach. The three loop 
BFKL predictions by three of the authors [13j [14l [HI [16] were used among other things in writing 
the symbol of the remainder function in general kinematics by Dixon, Henn and Drummond [17) . 
Their findings are in agreement with next-to-leading three loop calculations of one of the authors 
in collaboration with V. Fadin [18] . 

The Regge factorization made it possible to write the all loop LLA remainder function as a closed 
integral in Ref. [3] and calculate it analytically to all loops in a more restrictive kinematics (multi- 
Regge and collinear) by three of the authors in Ref. [Jj5] . This was shown to be consistent with the 
OPE expansion formula for the remainder function in the collinear limit derived by Alday, Gaiotto, 
Maldacena, Sever and Vieira [20] . An intriguing similarity of the OPE formulae (2jJ [22l [23l [24] 
and the BFKL expressions suggest some deep relation between these two approaches that describe 
very different kinematic regimes. These relation seems to be analogous to the relation between 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [25] and BFKL equations. 

The extensive use of symbol in the recent studies of the MHV and N fe MHV amplitudes made it 
possible to derive important analytic results in a series of very interesting publications [26 , 27j [28j 
[2"5] . The symbol for all leg two loop MHV amplitude derived by Caron-Huot in Ref. [3U] using an 
extended superspace is of particular interest to us because it is directly related to the main results 
of the present paper as follows. 

In the next section we consider the 2 — > 5 MHV scattering amplitude in a kinematic region, 
where Mandelstam cut gives a non-vanishing contribution. We call this region the Mandelstam 
region. Using the BFKL approach we derive the LLA remainder function in closed integral form 
of eq. ([28]) and calculate it analytically at two loops given by eq. ([36]). It turns out that the 
2 — > 5 remainder function can be written as a sum of two 2 — > 4 remainder functions due to 
the recursive properties of the leading order impact factors in the BFKL approach depicted in 
Fig. [T2] We generalize this result for the case of the 2 — > 2 + (n — 4) scattering and write the 
corresponding two loop LLA answer for an arbitrary n as a sum of six-particle remainder functions 
in eq. (j5Tj) and eq. (|56|) in a particular Mandelstam region. We choose this region such that we 
do not have contributions of Bartels-Kwiecinski-Praszalowicz (BKP) [30l EH] states that appear 
in other regions for n > 8. The Mandelstam regions with BKP states can be identified using 
dispersion representation of the 2 — ► 2 + (n — 4) similar to one in Ref. [2] • The number of terms in 
the dispersion representation is determined by the Steinmann relations 32 . 

The remainder functions in eq. ([36]) . eq. (f5T|) and eq. ([56]) are the main result of the present paper. 
The parallel study of one of the authors in collaboration with Spradlin, Vergu and Volovich [33] 
shows agreement between the BFKL calculations of the current paper and the symbol derived by 
Caron-Huot in Ref. |29] for the two loop remainder function with an arbitrary number of external 
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particles. 

Some relevant calculations are presented in the Appendices. 



2 2 — > 5 amplitude 

In this section we study the high energy behavior of the 2 — ¥ 5 scattering amplitude shown in 
Fig. [1] in Af = 4 supersymmetric Yang Mills (SYM) theory. In the leading logarithmic approxima- 
tion (LLA) only gluons contribute since in this limit i-channel exchanges are dominated by particles 
with highest spin. This fact allows us to apply the QCD-based BFKL approach to calculating the 
corrections to the BDS amplitude with the leading logarithmic accuracy. 



Pa Pa 




Figure 1: The 2 — > 5 amplitude. The produced particles ki are strongly ordered in rapidity in the 
multi-Regge kinematics. 

The Mandelstam invariants for the 2 — > 5 amplitude are given by 

s = (p A + p B ) 2 , s a > 12 3 = (pa> +h+k 2 + k 3 ) 2 , sb'123 = {PB> +ki+k 2 + k 3 ) 2 , (1) 
S123 = (ki +k 2 + k 3 ) 2 , sa'u = (pa> +h+ k 2 ) 2 , s B >23 = (pb> + k 2 + k 3 ) 2 , si = (pA> + fci) 2 , 

s 2 = (h + k 2 ) 2 , s 3 = (k 2 + fc 3 ) 2 , s 4 = (p B , + k 3 ) 2 , ti = (p A > - Pa) 2 , t 2 = (p A > + ki -p A ) 2 , 
h = (PA' +ki+k 2 - p A ) 2 , t 4 = (p A > + h + k 2 + k 3 - p A ) 2 = [pb 1 - Pb) 2 ■ 

The multi-Regge kinematics is characterized by a strong ordering in the rapidity of the produced 
particles 

S > SA'123, SB'123 > S123, S^'12, S B >23 > Si,S 2 ,S 3 , S 4 > -tl, -t 2 , -t 3 , -t A . (2) 

It is useful to introduce Sudakov parametrization 

fcf = aiV \ + (3 iPb + k? ± (3) 
for which the multi-Regge kinematics in eq. © reads 

1 > a x > a 2 > a 3 > 0. (4) 
The Pi components of fcf are also ordered due to the on-shellness condition 

-k 2 k 2 

ft = = ^, (5) 

OiiS CtiS 

where fc 2 ± = — k 2 . In the rest of the paper we use bold script to denote the transverse momenta. 
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The high energy scattering amplitudes have contributions of the Regge poles (poles in the 
complex angular momentum plane) and Mandelstam (Regge) cuts. In planar amplitudes the con- 
tributions of the Mandelstam cuts cancel out in most regions while in some regions (Mandelstam 
regions) they give logarithmically divergent terms (see Ref. [2 for more details). 

In this section we consider one of the Mandelstam regions for the 2 — > 5 scattering amplitude 
shown in Fig. [2] where we flip all three produced particles and have a non- vanishing Mandelstam 
cut contribution that corresponds to the discontinuity in s and Si23- In this region 

S,S2,S3,Sl23 > (6) 

and all other energy invariants are negative. 




Figure 2: Mandelstam region of the 2 — > 5 amplitude with a non- vanishing contribution of the 
discontinuities in s and S123 to the remainder function. 



For the scattering amplitudes with n external particles being on-shell we have 3n 
dent cross ratios 

~>2 ^,2 

2 2 
X kj X il 

written in terms of the dual coordinates defined by 

Pi=Xi-x i+ i. (8) 

In the case of the 2—^5 amplitude we have n — 7 and thus 6 independent cross ratios. However, it 
is beneficial for our purposes to follow the notation of the BDS paper and use 7 cross ratio, when 
we parameterize them in terms of the physical momenta. In the multi-Regge kinematics eq. Q the 



— 15 indepen- 



(7) 
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cross ratios we use are given by 



W73 = -2 — — = 7 77 7 - i) (y) 

^14^73 (-Sl23j(-SA'12) 
Xj]^ _ (-S123B')(~S3) _ 1 
^14^25 (-Sl23)(-S23B'J 

2, 



^62^35 = (Si)(-t 2 ) _ q2«3 _ ( 
^25^36 (-S23B')(-*3) ~ q 3 «2 ~ 

^73^46 = (s A 'i2)(-U) _ qik|a 3 

^36^47 (-SA'123)(-*3) ~~ q 3 k§a 2 



W47 - .0 -o = 7 — : 77 — : 7 - 1 



(-si 23 )(-s) 



(~Sl23S')(~ s ^'123) 



„ _ x 25 x l\ _ (-h)(-S23B>) _ qfog _ „ 
51 22 / j- \! \ — 2 — ' 

^51^62 (-*2)(-Sl23B') q 2 "l 

= ^63^72 _ (~Sl)(-*3) _ qjkl^g 
62 ^62^73 (-SA'12)(-t 2 ) q^k^ai 

In this kinematics we consider a symmetric subregion, where we have a 2 /ai ~ 0:3/0:2 = S for <5 
being a small parameter. Note that three of them go to unity with different rates 

(k 2 +k 3 ) 2 (k 2 + k 3 ) 2 a 3 

1 — 7il4 — — To 0, (10) 

S3 kg a 2 

(ki+k 2 ) 2 (k 1 +k 2 ) 2 a 2 
1 - m 73 ~ ~ -5 cx d, 

s 2 ai 

(k!+k 2 +k 3 ) 2 (ki + k 2 + k 3 ) 2 a 3 2 

1 — M47 ~ ~ — = CX . 

S123 kg o\ 

In the Mandelstam region of eq. ([6]) shown in Fig. [2] only one cross ratio has a non-zero phase 

u 47 = M e" i2 \ (11) 

It is worth emphasizing that this is also the only cross ratio that depends on the transverse mass 
(ki + k 2 + k 3 ) 2 of the bunch of the produced particles that are flipped in this Mandelstam region. 
It is also the most rapidly approaching the unity cross ratio as one can see from eq. I|10p. We use 
these two facts later to identify the cross ratio that has the phase for similar Mandelstam regions 
in the scattering amplitudes with a larger number of external gluons. 

Now we apply the BFKL approach to calculating the discontinuity in si 23 of the 2 — > 5 scat- 
tering amplitude to the leading logarithmic accuracy. In a way analogous to the 2 — !> 4 amplitude 
considered in Ref. [3] we can take advantage of the Regge factorization, which allows us to de- 
compose the amplitude into several blocks that we calculate separately. These blocks are shown 
in Fig. [3] The impact factors are the same as in the case of the 2 — > 4 amplitude. The BFKL 
Green function Gbfkl is universal and corresponds to the homogeneous octet-BFKL equation for 
the wave function / with removed propagators 

Ef = Hf. (12) 

The detailed discussion of this equation is presented in Ref. 3 and here we only need its eigenvalues 

and eigenfunctions 

/ u \ iw +§ / If* \ lv+ % 

wm) = (A) l-A^) ( i4 ) 
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with their completeness condition 

E / + °° duf^k 1 , q')UAk, q) = 2n 2 6 2 (k' - k) ^ ~ ^ . (15) 




Disc 



Figure 3: The 2 — > 5 amplitude factorized into blocks. <3?j represent impact factors, Gbfkl illustrates 
propagation of the BFKL state and C denotes the central emission block. The dashed line stands for 
the discontinuity in S123. 



The impact factor $ 2 in Fig. [4] 



$0 



k 3 (k"-q 2 ) 

q 3 {k"-k 2 -k 3 ) 



(16) 



was calculated in the terms of the anomalous dimension v and the conformal spin n in Ref. [3] 
(note a redefinition of momenta) 



X2 



2UI 



d 2 k 



2tt \k" - k 2 \ 2 \q 2 - k'f 



92 



1 1 



2i/i-f 

In a similar way one finds 



9 1 



ifi — m/2 



k" - k 2 , 
k 3 



qi 



k"* 



k 3 {k h 



92 J 



q 3 (k"-k 2 -k 3 ) 



(17) 



1 1 



Xi 



2 iv+ § 



9i 



-iv—n/2 



ql\- W+n/2 _ 1 (-1)" 
'fej 



— iv—n/2 



-iv-\-n/2 



.(18) 



A new piece in the 2 — > 5 amplitude compared to the 2 — > 4 amplitude is the central emission 
block C in Fig. EI In order to find the central emission block C we insert the completeness condition 
of the BFKL eigenfunctions between each two adjacent produced particles k{ as shown in Fig. [5] 
Each such insertion is denoted by a dashed line that also denotes insertions of the BFKL eigenvalue 
of eq. (Tj"3|) when we go to higher loops. 
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is* 



k"-k„ 



Figure 4: The impact factor <3>2 i n Fig- El 

























*** 












r-k 2 


k 3 



v, n 



Figure 5: Decomposition of 2 — > 5 amplitude into two impact factors and a central emission block. 
The dashed lines denote insertions of the BFKL eigenvalue. 



In the Appendix [5] we calculate the central emission block C illustrated in Fig. [6] in terms of 
the conformal spins n and m as well as the anomalous dimensions v and and the result reads 

(" =--( — ! I - I I - I ( ^ ) B, (19) 



2 \ ^2 



where / is given by 

b = (-iy 



F(l -iv- n/2) T(+ifi + m/2) T(i(v - u) + (m-n) /2) 



T(l +iv- n/2) T{-iu + m/2) T(l - i{v - u) + (m - n)/2) ' ^ 

It is worth emphasizing that C in eq. p^|) reduces to the form of the impact factors in eq. (fT7j) 
and eq. (fT8f in the following limits 



1 1 / „* \ ,Iy -™/ 2 / _ \ »f+n/2 

lim C m=0 = --7 I tt 

M->0 2 11/ - § Vfc| 



(21) 



f-^o 2 + § \& 2 



<?2 \ / 92 



After calculating all building blocks for the discontinuity in S123 of the 2 — > 5 scattering ampli- 
tude in Fig. [3] we convolve impact factors $1, $2 and the central emission block C together with 
the BFKL eigenvalues E Ufn and and obtain 



Reg Sl 



+00 +00 

J; ^ ; ^ / 

n— — 00 m=— co 



r+oo / \ -a(B„, n +E M , m ) 



(22) 
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k-k„ 



v, n 



|X, m 



Figure 6: The central emission block C in Fig. [5j 



The expression in eq. (|22l) is divergent at one loop and should be regularized. We fix regularization 
in Appendix [C] by comparing eq. (|22p to BDS amplitude at one loop. At higher loops this integral 
is finite due to the fact that E Vt „ — > for n = and v — > and all infrared divergencies are included 
in the gluon Regge trajectories u}(U) given by 

io(t) = a(^-ln { —pj, (23) 

where d = 4 — 2e and 

a s N c 



a 



2n 



(Ane-^Y . (24) 



Then the Mandelstam cut correction to the BDS amplitude at the leading logarithmic accuracy 
reads 

/ 04n + °° + °° r + oc r + ' x • ( /a \-*(.E„, n +En, m ) \\ 



n— — oo m— — oo 



The minus unity in the brackets on the right hand side of eq. (1251) removes the divergent one loop 
contribution. The choice of the energy scale so does not change the result in the leading logarithmic 
approximation and thus it is dictated only by the requirement of the dual conformal invariance and 
the Regge factorization. Among other possible choices we prefer 

£123 _ 1 gi |qa| [qgj |k 3 | ^ 

SO y/U25 "36 U51 «62 Qf 3 |qi | |q 4 | |ki | ' 

where the cross ratios Uij are given in eq. ([5]). So that the remainder function R defined by 

M = M BDS R 2 ^ 5 , R 2 ^ 5 = 1 + a 2 R { 2 % + a 3 R { 2 % + ... (27) 

reads 

where the product xi C xi ( see e Q s - (fl7l) - (fT9)) ) can be written as 

.. c .. (-i)" +m r(-w-|) r ( z M + f) r^- M) + ^) il/+f 

Here we introduce a complex variables io, expressed in terms of the complex transverse momenta 

kiq 3 k 2 q 4 , , 

toi = — — , w 2 = — ;— . (30) 
qik 2 q 2 k 3 

The variables Wi are cross ratios in the dual space of the transverse momentum. To this point we 
considered only the discontinuity in Si23- However there is also discontinuity in s for the 2 — > 5 
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Disc 

s 

Figure 7: Decomposition of the 2 — > 5 amplitude for the discontinuity in s. The amplitude is built 
of two impact factors $j, central emission block C and the BFKL Green function Gbfkl- 
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Figure 8: The lower impact factor 3>2 in Fig. 



scattering amplitude in the Mandelstam region shown in Fig. [2] The discontinuity in s of the 2 — > 5 
scattering amplitude is illustrated in Fig. [7J 

The only difference between calculation of the discontinuity in s and the discontinuity in S123 
is that the impact factors $i should be replaced by a slightly different impact factors For 
example, for i = 2 the impact factor for the discontinuity in s is depicted in Fig. [8] 

The impact factor $2 was calculated in Ref. [3] 



qi(k" - k 2 ) 



Due to the fact that 



1 



(31) 



(32) 



the discontinuity in s is expressed through integral we had before for the discontinuity in S123, 
though with a different regularization at one loop in full analog with the 2—^4 case. Thus starting 
at two loops the contribution of the discontinuity in s to the remainder function is the same as of 
the discontinuity in S123. 
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The LLA remainder function i?2->5 in eq. (|28|) at two loops to the leading logarithmic accuracy 
was calculated in Appendix [Cl 



#2_U = y In V u 25 "36 U51 "62 I In 1 1 + w 2 + WXW2 1 2 In 



1 + W 2 + Wl»2 



w 2 (l + Wl) 



(33) 



In 



I + W2 + W1W2 



W1W2 



In 



1 + w 2 + W1W2 



1 + W2 



«7T / S123 

— In 



+ ln 



ln qg(ki +k 2 + k 3 ) 2 ^ q 2 (ki + k 2 + k 3 ) 2 



n 2 k 2 
Vll K 3 



ql(ki+k 2 ) 2 



q!(ki + k 2 + k 3 ) 2 ^ _q 2 (ki + k 2 + k 3 ) 2 



q 2 (k 2 +k 3 ) 2 



The remainder function i? 2 4-5 nas an interesting property of being expressed through the cor- 
responding remainder function for the 2 — » 4 amplitude found in Ref. [3] 



/■v., >• 



(34) 



where 



/e(w, u>*) = In J 1 + w\ 2 ln 



1 

1 + - 

w 



(35) 



We can compactly write the expression in eq. ()33|) as a sum 



#2M-5 = y InV^ «36 «51 "62 (fe{w a ,wl) + f 6 (w b ,W%)) , 



(36) 



where 



«'<i 



1 1 



">b = 

1 + — W 2 1 + W'l 



(37) 



In the next section we generalize the present discussion to the 2 — s- 2 + [n — 4) scattering 
amplitudes with an arbitrary number of produced gluons n — 4. 



3 Some generalizations for more legs 

In this section we apply the results of the previous section to the 2 — > 2 + (n — 4) scattering 
amplitudes with n — 4 produced gluons illustrated in Fig. |H1 We argue that at two loops the 
corresponding remainder function in the Mandelstam regions, where we flip at least two adjacent 
particles can be written as a linear combination of the 2 — ¥ 4 remainder function. This happens 
due to a very special structure of the effective emission for a definite helicity vertex of the produced 
particles 

V2^f (38) 
k 

shown in Fig. [TOl 

If we consider an effective diagram with emission of two gluons with the same helicity, when 
the transverse propagator between the two emissions is canceled as follows 

^2 ^^^2^ = 2^4. (39) 
ki I92I 2 &2 ^1^2 
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n-3 



n-4 



Figure 9: The 2 — > 2 + (n — 4) scattering amplitude with n — 4 produced gluons. 



Figure 10: The effective reggeon-particle-reggeon emission vertex. The t-channel gluons are not yet 
dressed at the Born approximation and thus are denoted by simple (not curvy) lines. 



Using this property we can find the impact factor $34 with two emissions for the discontinuity 
in s, similar to $2 of eq. (|3"TT) . Let us recall how one calculates impact factor $2 depicted in Fig. [5] 
Plugging in the effective emission vertex and the transverse propagator we write 



V2 



(k" - k 2 )(k" - k 2 - k 3 y 



1 



k* 3 \k"-k 2 -k 3 \ 



and then divide this by Born expression 
-(<72 - k 2 )(q 2 -k 2 - k 3 )* 



|<?2 - k 2 - k 3 



v/2 



y/2- 



k* 3 {k" - k 2 - hi) 
(<72 - k 2 ) 



As a result we have 



$2 



('12 



h)(k" 



k 3 (q 2 - k 2 - k 3 ) 
qA(k"-k 2 ) 



(<?2 - k 2 ){k" -k 2 - k 3 ) q 3 {k" -k 2 - k 3 ) ■ 



(40) 



(41) 



(42) 



Next we introduce one more produced gluon k^ as illustrated in Fig. [TT] and calculate the 
corresponding expression by plugging in the effective vertices and the transverse gluon propagators 



V2 



(k" -k 2 )(k" -k 2 -k 3 y 



^3 

= 2- 



k 3 \ 2 



V2 



{k" -k 2 -k 3 )(k" -k 2 -k 3 -ki) 



k^ 



fc 3 -fc 4 | 2 



1 

k^. kr> k 



(43) 



3^4 
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k"-k„ 



k"-k 2 -k 3 -k 4 



Figure 11: Impact factor for the discontinuity in s with two produced gluons. 



Dividing this by the corresponding Born expression 



V2 



(<72 - k 2 )(q 2 -k 2 - k 3 )* 



1 



(12 



\q 2 - k 2 - k 3 
1 



rV2 



(q 2 -k 2 - k 3 )(q 2 - k 2 ~ k 3 - fc 4 ) 



\q 2 - k 2 ~ k 3 



q 2 — k 2 — k 3 — ki k 3 kl 



we readily get 



$34 = 



(g2 - k 2 - k 3 - kj){k" - fcg) = 

(q 2 - k 2 )(k" -k 2 -k 3 - fc 4 ) q 3 (k" - k 2 - k 3 - k 4 ) ' 



q 5 (k"-k 2 ) 



(45) 



An important observation is in order. The impact factor for an emission of two gluons with the same 
helicity can be obtained from a corresponding impact factor with one gluon emission by shifting 
the transverse momentum of the produced gluon. In other words impact factor $34 is the same as 
$2 with k 3 —¥ k 3 + k±. Graphically this statement means that any two adjacent gluon emissions 
of the same helicity can be represented as a stretched diagram in the transverse space, where the 
two gluons are emitted at the same point as shown in Fig. 1121 In a more general case of many 



k"-k„ 



*2# 




k"-k 2 


















k"-k 2 -k 3 -k 4 









«4r 

(44) 



Figure 12: The impact factor $34 for two produced gluons of the same helicity with momenta k% and 
&4 can be written as an impact factor with one produced gluon with momentum k^ + k^. 



gluon emissions of the same helicity all of them effectively emitted from the same point due to the 
cancelation of the propagators between them as shown in Fig. 1131 The corresponding impact factor 
$ m for m emitted gluons can be obtained from $2 in eq. pip by shifting the transverse momentum 
&3 — > &4 + ••• + k m . This fact allows us to factorize any 2 — > 2 + (n — 4) amplitude into two pieces 
with redefined impact factors. These two new redefined impact factors are then convolved with the 
BFKL propagators in a way it was done for the 2 — > 4 amplitude and thus give the 2—^4 answer 
with redefined momenta. The number of factorization points determines the number of 2 — > 4-like 
terms in the final answer. For Mandelstam region of the 2 — > 2 + (n — 4) scattering amplitude, 
where all n — 4 produced particles are flipped we have n — 5 factorization points denoted by a 
dashed line in Fig. 1141 In the case of the 2—^5 amplitude we have two factorization points, which 
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Figure 13: The impact factor for % + ...+ j produced gluons of the same helicity with momenta k%, 
fci+i, ••• ,kj can be written as an impact factor with one produced gluon with momentum fcj + ... + kj. 

explains why the two-loop remainder function in eq. (|36p can be written as a linear combination of 
two 2 — > 4 remainder function. 




Figure 14: The 2 — > 2 + (n — 4) scattering amplitude can be factorized into two pieces, which are 
then summed over. The dashed line denotes the factorization point. 

It is worth emphasizing that here we considered impact factors $i for the discontinuity in s. 
We showed (see eq. (13"2")) and text wherein) that in the case of the 2 — > 5 amplitude they differ from 
those for discontinuity in S123 only by an additive number 

$i + $i = l. (46) 

The difference between them is important only at one loop level and changes the regularization 
prescription. Starting at two loops it vanishes after the integration over the loop momenta because 
of the BFKL propagators. The best way to see this is to look at the integral representation eq. (|22"1) . 
which is divergent at one loop at v = (/J, = 0) for n = (m = 0), while starting at two loops this 
divergency is absent due to the fact that the BFKL eigenfunctions E v _ n vanish at those points. Any 
emission of one or more produced gluons does not change this property and the self consistency 
requires the remainder function for the 2 — > 2 + (n — 4) MHV amplitude in this Mandelstam region 
calculated from the discontinuities in Si...(„_4) and s to be the same. 

This simple iterative structure breaks down if we go beyond two loops. However, we expect a 
similar, though more complicated iterations at higher loops, which will be published by us elsewhere. 

At two loops we can formulate a simple mnemonic rule based on the iterative property of the 
impact factors outlined above. For any given 2 — >• 2 + (n — 4) scattering MHV amplitude in the 
Mandelstam region, where all n — 4 produced particles are flipped, the amplitude can be factorized 
into two pieces in n — 5 ways as shown in Fig. Q3] The dashed line denotes the factorization point. 
Let us consider the upper piece in Fig. [14] in more details. We insert the completeness condition 
eq. (Ti"5|) for the BFKL eigenfunction in the factorization point and integrate over the transverse 
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momenta in a similar way we did in eq. (|17|) . The impact factor $123 in v, n space for the emission 
of three produced gluons with momenta k\ , and k^ depends only on the following ratio of the 
complex transverse momenta 

qi -ir- (47) 



ki + k 2 + k 3 

In a similar way for the lower piece in Fig. [T?]for n — 7 produced particles with momenta fc 4 , 
fc„_5 and k n _4 one has (note inverse momenta ordering) 

ki + h + ■■■ + k n -4 , AO , 

Qn-3 

The product of two impact factors depends only on the complex cross ratio in the space of the 
transverse momenta 

gl (fc 4 + fc 5 + ... + fc»_ 4 ) 
3 ~ q n -3(k 1 + k 2 + k 3 ) ' l4yJ 

where the subindex of W3 stands for the factorization point between gluons with momenta k% and 
fc 4 . In general, it is given by 

Mt = ft(^i+ : -+ fc ;-*> . (50) 

q n -3(ki + ... + h) 

Then the two loop remainder function for the 2 — > 2+(n— 4) scattering amplitude in the Mandelstam 
region, where we flip all of the n — 4 produced particles reads 



n— 5 



The function of the transverse momenta fe(w,w*) appears in the remainder function for 2 — > 4 
scattering amplitude and was defined in eq. (|35l) . Using the transverse momenta conservation 

qi = q n -i + fci + ... + fc„_4 (52) 

we write (see Appendix [D| 
f( ^ 1 |3i 1 , 1 2 1 q?+i( k i + - + k "-4) 2 1 q 2 +1 (k 1 + -+k„_ 4 ) 2 

/ 6 K^) = in|i + ^| mi + - qL3(kl + ,.. + kt)2 % 2(ki+1 + ... + kn _ 4)2 - ( 53 ) 

In the Mandelstam region under consideration, where we flip all produced particles, only one 
cross ratio processes the phase — i2ir, namely the one that has both s and si...„_ 4 in the numerator. 
It is easy to identify it because it is goes to unity with the fastest rate in the multi-Regge kinematics 

l-U~5 n - 5 , (54) 

where 

<5=^±i (55) 
m 

is defined in terms of the Sudakov variables a.i in eq. (J3J . 

In the next section we discuss other Mandelstam regions of the 2 — » 2 + [n — 4) scattering 
amplitude at two loops. 
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3.1 Other Mandelstam regions at two loops in LLA 



We can also consider Mandelstam regions of the 2 — > 2 + (n — 4) scattering amplitude, where we 
flip less than n — 4 adjacent produced gluons having the same helicity. It is worth emphasizing 
that this discussion is not limited to MHV amplitudes and the only condition is that the produced 
gluons we flip must have the same helicity. 

Any gluon emission of any helicity that does not participate in building the Mandelstam cut is 
factorized out and canceled when we divide the impact factor for the corresponding Born expression. 
To illustrate this statement let us consider 2 — » 6 scattering amplitude in Fig. [5] in the Mandelstam 
region where we flip first three produced particles with momenta k\, k 2 and ^3. In this region we 
find the discontinuity in S123 and, using the notation of the previous section, we can have only two 
factorization points between k\ and fc 2 , and between fc 2 and fc 3 . The factorization point between 
&3 and A)4 is excluded because it does not correspond to the discontinuity in S123. 

Let us look closely at the lower impact factor for the factorization point between gluons with 
momenta k 2 and k^ shown in Fig. 1151 Due to the Regge factorization any dependence on k^ and q§ 



q^ 



k"-k n 



Figure 15: The impact factor with emission of two gluons with momenta and £4. The helicity of 
the produced gluons is not necessarily the same. 



in the impact factor in Fig. [15] cancels if we divide it by the corresponding expression in the Born 
approximation depicted in Fig. 1161 As a result we get back the same expression we had for $2 in 
Fig.H 

This means that to the leading order the impact factor that stands for Fig. [T5l scaled by its Born 
expression is exactly the same as $2 that corresponds to Fig. [4] This statement is independent on 
the helicity of the produced gluon with momentum k^ and valid for any number of gluons of any 
helicity produced below the one with momentum k%. The same is true if we add gluon emissions 
above the first gluon that spans the energy for which we calculate the discontinuity of the amplitude. 

Thus we can formulate a leading logarithmic prediction for the remainder function for an arbi- 
trary MHV 2 — >• 2 + (n — 4) amplitude at two loops in the Mandelstam region, where we flip any 
k — to adjacent produced particle 

(\ ™} 
-^J^2f 6 {wi,w*), k>l, m<n-4. (56) 

The indices k and m label the first and the last particles that span the energy s&.. m for which 
we calculate the discontinuity. The only condition here is that all of the flipped particles should 
have the same helicity, while the helicity configuration of all other produced particles is arbitrary 
and does not effect the remainder function in eq. (|56]l. The phases of the cross ratios for these 
Mandelstam regions are calculated case by case. 
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k 3 



Figure 16: The Born expression for the two gluon production. The t-channel gluons are not yet 
reggeized in the Born approximation and the curvy lines are introduced merely to distinguish between 
reggeons and particles. 

All results of this section are valid only in the leading logarithmic approximation. The next- 
to-leading corrections to the effective gluon emission vertex break the recursive properties of the 
impact factors. 

4 Conclusions 

We calculated the remainder function for the 2 — > 5 scattering MHV amplitude with leading 
logarithmic accuracy in the Mandelstam region where s, S2, S3, S123 > that corresponds to all three 
produced particles being flipped as shown in Fig. [2] The result is given in the integral representation 
of eq. (|28|) to any loop order and calculated analytically at two loops (see eq. (|36|)). We found that 
the two loop leading log remainder function in this Mandelstam region can be compactly written 
as a sum of two remainder functions for the 2 — > 4 scattering amplitude calculated in Ref. 3 . This 
iterative structure happens due to the fact that the impact factor for two or more produced gluons 
is expressed though the impact factor for one gluon emission with a shifted transverse momentum. 
This property of the impact factors makes it possible to derive the remainder function for the 
2 —> 2 + (n — 4) MHV amplitude in the Mandelstam region, where all n — 4 produced gluons are 
flipped and have the same helicity (see eq. ([FTjO . In this region we do not have contributions of the 
Bartels-Kwiecinski-Praszalowicz (BKP) states, which appear in other regions for n > 8. In that 
sense the computations are the similar to the n = 6 case. This explains why the 2 — > 2 + (n — 4) 
result can be written as a sum of n — 5 LLA remainder functions for n = 6 MHV amplitude. 

Furthermore we consider other Mandelstam regions of the 2 — >• 2 + (n — 4) amplitude, where we 
flip any number of adjacent produced gluons having the same helicity. The corresponding remainder 
functions at two loops are given by eq. (f56|) . The last result is also valid for N fc MHV amplitudes, 
where the nipped particles have the same helicity, while the helicity of all other produced particles 
is arbitrary and does not effect the contribution of the Mandelstam cuts to the leading logarithmic 
accuracy. 
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A Central emission block 



In this section we calculate the central emission block shown in Fig. [5] The dashed lines denote 
the insertion of the completeness condition for the BFKL eigenfunctions eq. (fT5j) and the emission 
of the gluon with definite helicity is given by the effective emission vertex in eq. psp . Thus we can 
write for the central emission block in Fig. [5] the following integral 



J 



d 2 k 



-iv-\-n/2 



k 



Q2 — k 



-iv — n/2 



k{^k ) 



\k\ 2 \k - k 2 \ 2 \q 2 - k\ 2 k* \&-k 



i[i—m/2 



k - k 2 
92 - k 



z/i+m/2 



It is useful to introduce the dual coordinates in the transverse momentum space to exploit the 
conformal properties of J. The dual coordinates are depicted in Fig. [TjJ and defined by 



92 - k 2 = z B - zo, q 2 = za - z , k 2 = z A - z B , k = z A - z > 
k - k 2 — z B - Zo> q 2 - k = zq> - z . 



(A.2) 



(A.l) 




Figure 17: Dual coordinates for the integral representation of the central emission block. 



The behavior of the integral J under inversion allows us to factor out the dependence on the 
external momenta 



(Pzp 

IT 

inversion 



'- {zAvT iv+ ^ x M~ iv ~^ &r _f (W** -1 04 o r~ f ~ iM+f _1 (^' ) i!/+f - iM " 



d 2 z 



L {ZAvT iV+ ^ 1 {ZAvT iV ^ (^0')^ f (ZBO'Y^- 1 (Z* , ) 



iv—^—ifjL+^—l 



(zo'o) 



—i^i- 



Thus we can write 

1 q 2 (q* 2 -k*) 1 



J 



ml 



k 2 \q 2 ~k 2 \ 2 \q 2 -k* 



92 - h 



92 



The factor B is a c-number, which is obtained by taking q 2 — > oo and k 2 = 1 in J 

f d 2 k 

B = n„*\—iv+n/2-l fe-iu-n/2 _ ^yfj,-m/2 ^ _ ^yn+m/2-1 



To calculate B we introduce the Sudakov variables 

k = k\ + ik 2 — x, k* = k\ — ik 2 = y 



B (A.4) 



(A.5) 



(A.6) 



(A.3) 
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and write 



B= f dxdy x~ n+1 (x - I)'"' 1 

J 2-ki (xy - ie) i "- 2 5+ 1 ((1 - x)(l - y) - ie)"V+? 1 

1 i7r(— 1) _ p i7r(— — 1) 

dra - *""'^ - aO^^C-l)" 1 " 1 

27TI 

_! T(l - if - n/2) r(+i/i + m/2) r(i(i/ - /i) + (m - n) /2) 



(-l) r 



T(l + ii/ - n/2) r(-ip + m/2) T(l - i(> - fjt) + (to - n)/2) 



Then the central emission block C is obtained by dividing J by the corresponding Born expres- 
sion 



1 g 2 (g 2 *-fc 2 *) 1 
|g 2 | 2 k* \q 2 -k 2 \* 

and reads 

1 / 93 A / 93 \ 2 / 92 \ 2 / 92 



(A.8) 



C =-AW VI) [fJ ( -> B (A - 9) 

We normalize C in a way consistent with limits in eq. (|21[) . 



B One loop from BDS in the Mandelstam region 

In this section we consider the BDS formula at one loop for fixing the normalization of the remainder 
function i? 2 ->5 ■ We perform the analytic continuation of the one loop expression to the Mandelstam 
region shown in Fig. [2j where we flip all produced particles with momenta fcj. In this Mandelstam 
region we have s, s 2 , S3, si 2 3 > and only one cross ratio that has a phase is 

u 47 = \u 47 \ e~ a \ (B.l) 

Firstly we decompose the BDS formula into two pieces - one corresponding to Regge poles and 
the other coming from the beginning of the Mandelstam cut at one loop. Generally, the piece that 
stands for Regge poles for 2 — > 2 + (n — 4) amplitude can be written as follows 

™2 ) r_Rp fi (t„_ 4 ,i„_3,in(-K n _4)) f n 2 ) rp_ap(t n _ 3 ), (B.2) 

where the particle- reggeon-particle Tp RR and reggeon-particle-reggeon T R p R vertices were found 
in Ref. [5] from the BDS amplitude with n = 4 and n — 5 external gluons. To the first order in the 
coupling constant they are given by 

lnr PH p(t) = -i + i-ln^+2C2 (B.3) 
t l 2e /ir 



and 

1 1 9 — k 1 — k / (-ti)(-< 2 ) 1\ 1 , o -*i 1, 
IB W*i, fe, In -«) = - ^ - i in 2 ^ + - In (in _ _ j _ _ _i _ _ & (B . 



The gluon Regge trajectory is 
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The parameters Ki depend only on transverse momenta in the multi-Regge kinematics and for the 
2 — > 5 amplitude are given by 



Kl = — ; — , -K 2 = —, r— , -«3 = -7 (B.6) 

\-SA'Yl) (-S123) (-S23B') 



Then we subtract from the n = 7 BDS amplitude the Regge pole contribution eq. (|B.2j) at 
one loop and get 

P 1 1 t. 1 1 1 ( - 32 ^ 2 L 1 1 1 L 1 1 1 (~ S 123)M 2 

i* = — In U47U73U14 + — hi U73 In — ; + — In uu In -p — ; + — In U47 in 



2e ~ ™* 1 2— "-(-txX-ta) ' 2"""™(-t 2 )(-t 4 ) ' 2 (-i 4 )(-ti) 
7r 2 In 2 Ux4 In 2 U47 In 2 M73 1 



3 4 4 4 2 



5^Li 2 (l-« i , i+3 ). (B.7) 



1=1 



The function i* 1 carries the information about the Mandelstam cuts and is given in an arbitrary 
kinematics. This fact allows us to perform the analytic continuation of F to the Mandelstam region 
where U47 = IU47I e _l27r in the multi-Regge kinematics and obtain 

P »tt . , (ki+ka+ka) 2 ^ 2 

i< ~ mm 5—5 . (B.8) 

e QIQ4 

Note that F — > for U14, U73, M47 — > 1 in the multi-Regge kinematics. 

C Integral at one and two loops 

In this section we calculate the integral in eq. (|2"51) 

+00 +00 ^.+00 r+00 
lW = E E / dv dn xi C X 2 (E„, n + E^y- 1 (C.l) 

for I = 1 and ^ = 2. 

The integrand is given by 

XlCx2 = 8 r(i + - f ) r(i ^Tf) r(i - ifr - M ) + ^) Wl K) " 2 K) 

where is given by eq. (fT3|) 

= -^^^ + ^ + iv + + il> (l - ii/ + - 2^(1) (C.2) 
and u>i are expressed in terms of the complex transverse momenta 

wi = — j~, W2 = — r~. (C.3) 
It is useful to introduce the phase and the square of the absolute value of Wi as follows 

a t = ^, fr = H 2 . (C.4) 

Firstly we calculate I £ at one loop I = 1 

7 = E £ ^ 8 r ( i + ^-f)r ( i- v + f ) (c - 5) 

n=— oo m=— oo ^ ^ v 2 /vr 2 / 

r(t(l/ - /i) + ^) n/Uf*, m/9fAn 

X T{l- l ( V -n) + ^f 1 Pl " 2 ^ 2 ' 
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The calculation of the integral in eq. (|C.5|) becomes much simpler if the consider only one region, 
where we have /3i < 1 and /?2 > 1 which is consistent with wi -H- l/w2 symmetry. Moreover we 
focus on summing contributions of yftijcti in that region and then exploit the Wi f-> w* symmetry 
to restore the full answer. 

The three gamma functions in the numerator of I^ 1 ' have an infinite number of the simple poles. 
We show that only first two gamma functions contribute in the region under consideration and that 
out of an infinite number of poles only one has non-vanishing residue. This is consistent with the 
2 — >■ 4 case, where integrand that corresponds to one loop expression has only one simple pole. 

The poles of the first gamma function in the numerator are located at 



n 

lV ~2 



a = 0,1, 



and the second gamma function has poles at 

m 

Thus the contribution of the residues at these poles is given by 



-t, t = 0,l,.. 



-l 

E 

n=—oo 771— 1 



E 



(-l) n+m T{m-n + s + t) 



r(i 



s )r(i 



t)T(l 



t) 



a 



/2 R -m/2 
2 l J 2 i 



(C.6) 



(C.7) 



(C.8) 



where the limits of the summation arc determined by the convergency of the integrals over v and 
/i at the large circle in the complex planes for fii < 1 and f3 2 > 1- There are also contributions at 
n = and m = 0, which we consider separately. We see that T(l — s — t) makes this expression to 
vanish for all s and t except s = t = and therefore it becomes 



7T" 

T 



E E 



(-!)"+"' T(m-n) n/2 ^ n/2 m /2 




0: 



-m/2 



In 1 



(C.9) 



Next we find the corresponding terms for n = 0, m ^ 



2 +°° 
7T x - 



r(m) m/ 2 R - m /2_ 
<°2 



2 ^ r(i 

™-1 v 



-2- In 
2 



1 + 4/^ 



02_\ 



(C.IO) 



and for m = 0, n 7^ 



4t 



HTTH n/2 R -n/2_ 



r(i 



Pi 




Adding these two to eq. (|C.9|) we obtain 



(C.ll) 



(C.12) 



Going back to the integrand in eq. (|C2[) we see that we have also poles from the third gamma 
function in the numerator, namely from T(i{v — /i) + T "~" ). These poles lead to a system of 
inequalities for the arguments of the gamma functions in the denominator of the integrand in 
eq. (|C.2I) that comes from the requirement of the arguments being greater than zero for non- 
vanishing integrand. Those inequalities have no solution in the region under consideration. Going 
back to the complex variables Wi we write 





•In 1 



1 

Wo 



(C.13) 
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and the symmetrization in Wi -f> w* gives 



^2 



(l) 7r " u.i.„.. . 1 - ^ . q^(k!+k 2 +k 3 ) 



r ; -__1ti 



1 

1 + tox H 

w 2 



= _:_ ln ^^"y "V J/ ■ (C.14) 
2 q 2 q 2 k 2 1 y 



For n — m = the integral diverges at /i = and need to be regularized. For (3i < 1 
and /?2 > 1 we choose a regularization that is compatible with the BDS amplitude as follows. By 
comparison with eq. IC.14I we get 

^oo ^ ^ rMrH^r^-i/.) tt 2 / ^ q 2 q 2 \ 

Reg -L dM s m- ¥ m + i,ni-iu + i,) = Y{-T +ln w^)' (CU5) 

which is in an agreement with the regularization for the 2 — ¥ 4 amplitude found in Ref. [3]. 
Next we calculate at two loops 



n— — oo m— — oo 



X r(l-i( I /- M ) + ^) ai ^ 2 ^ (^,n+^,n.). 

The two loop expression i' 2 ' is finite at n = m = and v, /i — due to the fact that the BFKL 
eigenvalues E v>n and E^^ m in eq. (|C.2[) vanish at those points. We also note that the symmetry 
n o —to, ^ <-> — ^ of I^ 2 -* corresponds to wi O l/w<i. Thus it is enough to find only the contribution 
of -Ev, n and then that of E^ }m in eq. (|C.16[) is obtained by symmetrizing the result with respect 
to w\ -H- l/w%. Moreover taking into account that the result should be symmetric in Wi it)*, 
it is enough to calculate only the piece that depends on w* — y/ pi/ai. As in the case of the 
one loop calculation we pick up the region that is compatible with w\ ■<-> 1 lw% symmetry, namely 
Pi < 1, fi% > 1. We calculate 

+oo + + +m rH ,_n) r(» M +f) 



n— — oo m— — oo 



r(i(f - M) + I V I ) n/2 w m /2 ^ 

" =^TT a l Pi a 2 "2 ^.n- 



r(i-^- M ) + ^) 



The integration over /i and the summation over m is not difficult to do if we recall that at one loop 
all poles that have non- vanishing residue in the our region come from T(i/i + tt), while all residues 
of poles of T(i(v — n)+ m ~ n ) in /1 are zero. The integrand of eq. (|C22|) differs from that of the one 
loop integral in eq. (|CJ.5[) only by E un , which does not effect the integration over ji. As we saw in 
the one loop case, the poles of T(i/j, + ^) are located at ifi + ^ = —t, t = 0,l, but due to the 
very special structure of the integrand only pole with t = contribute. Using the Cauchy theorem 
we integrate eq. (|C.22|) over fi and obtain 



-00 +00 



2n ^ ^ / dv - r(1 + w _ f) T{l _ w _n ) r(TT^) Ev ' n (c - 18) 



+00 ,.+00 \n n % oiv / n — \ 

dv { - I "'J 1 E v Jl + 



4 ^ 

n— — 00 

The expression in eq. (|C.18|) is reduced to the two loop integral for the n = 6 remainder function 
calculated in Ref. [3]( for j3\ < 1 and fi% > 1) 

J E f + °° d ^-^lt £,..„ : -br | I \/- j -1,, ! I \/- | In ; K '.!!)) 
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with redefined variables 



a = a\ I 1 



1 + 



(C.20) 



which correspond to 



Symmetrizing the result in W4 ■<-> «;* we get 



+00 +00 

E E 

n— — 00 m= — 00 



(lis 

r(i(i/ - m) + 



(-i) r 



r(^ + 



r(i + ti/-§) r(i-i/x + f ) 



(C.21) 



(C.22) 



T(l-i(v-fi) + 



-n \ 2 



2 

~j2 



W\ 



In 



1 + - 



The contribution of £^, m is obtained by substitution Wi -O- 1 /t/^ and we finally get the two loop 
answer for the n — 7 amplitude in terms of the function fe(w, w*) that appears in n = 6 amplitude 



^ = j(/6K,<) + / 6 K<)) 



where 



and 



/e(iy,ty*) = In |1 + iu| In 



1 



1 



1 

W2 



W 2 1 + Wi 



(C.23) 



(C.24) 



(C.25) 



D Larger number of external gluons 

In the previous section we found that the n = 7 remainder function can be expressed in terms of 
the n = 6 remainder function to the leading logarithmic accuracy in the Mandelstam region under 
consideration. The energy dependence is known and the problem is reduced to calculating the 
finite function of transverse momenta. In section |3] we showed that the leading order impact factors 
have a recursive properties for emissions of the gluons with the same helicity. This happens due 
to the effective cancelation of the transverse propagators between the emitted gluons, and thus the 
impact factor with m emitted gluons can be written as an impact factor for one gluon with shifted 
transverse momentum ki + fci+i + ... + k m . This property allows us to factorize 2 — > 2 + (n — 4) 
amplitude with n — 4 produced gluons in n — 5 ways. In the Mandelstam region, where we flip all 
of the produced gluons we have 



gi(fcj+i + ■■■ + kn-4) 
q n -3(h + ... + h) 



(D.l) 



where index i denotes the factorization point between gluons with momenta i and i + 1. Using the 
transverse momenta conservation 



k% + ... + fc„_ 4 - q n -3 



we can calculate 



1 + m = 



g n -3(h + ... + ki) + gi(k i+1 + ... + fc„_ 4 ) 
q n -3(ki + ■■■ + h) 



(D.2) 



(D.3) 
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The numerator of eq. (|D.3[) is simplified as follows 



q n -a(fa + ■■■ + fa) + qi{fa+i + ■■■ + K-i) = fe+i - fa+i - ■■■ - k n -4,)(fa + ... + fa) 

+ {q l+ i + fa + ... + fa)(fa + i + ... + fe„_ 4 ) = qi+iifa + .. + fc„_ 4 ) (D.4) 



and we get 



We also need 



1 + Wi = 



gj+i(fci + ■■ + k n -j) 
q n -z(k\ + ... + fa) 



1 J_ _ qi+i(fa + ■■ + k n -4.) 



(D.5) 



(D.6) 



Wi q\{fa+i + — + k n -i) 
and finally write the corresponding n — 6 remainder function in terms of the transverse momenta 



f 6 (wi,w*) = In |1 + i0i| 2 ln 



1 



f 



In 



qi+i( k i 



kn- 4 ) 2 q 2 + i(ki 



k„_ 4 ) 2 



q 2 _ 3 ( k i + - + k 4 ) 2 qi( k ^+i 



k „- 4 ) 5 



(D.7) 



This result is valid only for the 2 — » 2 + (n — 4) amplitude in the Mandelstam region, where we flip 
all 7i — 4 produced gluons. 
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